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ABSTRACT 

Aims. We provide a guideline to interpret the UVCS emission lines (in particular O VI and Si XII) during shock wave propagation in 
the outer solar corona. 

Methods. We use a numerical MHD model performing a set of simulations of shock waves generated in the corona and from the 
result we compute the plasma emission for the O VI and Si XII including the effects of NEI. We analyze the radiative and spectral 
properties of our model with the support of a detailed radiation model including Doppler dimming and an analytical model for shocks, 
and, finally, we synthesize the expected O VI 1032A line profile. 

Results. We explain several spectral features of the observations like the absence of discontinuities in the O VI emission during the 
shock passage, the brightening of Si XII emission and the width of the lines. We use our model also to give very simple and general 
predictions for the strength of the line wings due to the ions shock heating and on the line shape for Limb CMEs or Halo CMEs. 
Conclusions. The emission coming from post-shock region in the solar corona roughly agrees with the emission from a simple planar 
and adiabatic shock, but the effect of thermal conduction and the magnetic field may be important depending on the event parameters. 
Doppler dimming significantly influences the O VI emission while Si XII line brightens mainly because of the shock compression. 
Significant shock heating is responsible for the wide and faint component of the O VI line usually observed which may be taken as a 
shock signature in the solar corona. 
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1. Introduction 

Strong shock waves are expected to develop from fast Coronal 
Mass Ejections (CME). According to the most popular idea 
about the structure of CMEs, an expanding magnetic flux rope 
propagates in the solar corona and this leads to the creation of 
the CME leading edge and to the creation of shocks, as well. It is 
widely accepted that the front side halo CMEs are stroiig ly cor- 
related wit h geomagnetic disturbanc es (Brueckner et al. 1998; 

gOOO; .Gopalswamv et al.1 .2000.; Webb et al. 200C : 



Cane et al 



Zhang et al 



2003b . Thus, flie detection of CME related shocks 



is a crucial topic in space weather research. However, it is also 
a difficult problem, because it cannot be confirmed from white 
light images alone. A further problem is to describe the shock 
and figure out whether the environment or the CME core ex- 
pansion determine the shock wave properties, since they influ- 
ence the acceleration of energetic particles as the shock evolves, 
which is again a key question in space weather research. 

All proposed hypotheses agree that the shock wave, when- 
ever exists, stays on the leading edge of the CME (' HundhausenI 
[l987; Wagner & MacOueen 1983) though it may be too faint to 
be easily seen. The leading edge, as observed in white light im- 
ages, could be either plasma compressed by a shock or denser 
material in magnetic loops ejected from near the solar surface. 
Several shocks have been detected during CME front propa- 
gation with different diagnostics (Ciaravella et al. 2006). Sharp 
edges in white light images have been taken as evidence of 
shocks at the CME leading edge, but it is impossible to prove 



that they are shocks without the s upport of deeper diagnostics. 
For instance rVourUdas et al.l ( 1200 3) support the identification of 
shocks in SoHO/LASCO observations thanks to an MHD model 
where a shock was needed to fit the observational constraints. 
Sometimes the shock is detected without uncertainness by the 
presence of a Type II radio burst that gives a reliable estima- 
tion of the post-shock density, but not of speed and location. 
However, spectral analysis is the most general way to detect and 
describe the shock in the solar corona. 

The U l travio let Coronograph Spectrometer (UVCS) 
(iKohl et al.l Il995h o n board the Solar an d Heliospheric 
Observatory (SoHO) (Domingo & Poland 1988) gave a strong 
impulse to this research topic. Ciaravel la et al. (2006) pointed 
out the correlation between the shocks and several spectral 
properties of the post-shock emission. Overall, they found 
that the shock emission is comparable or weaker than the 
background coronal emission. The O VI 1032 A line is usually 
broad, but its intensity does not increase. Instead the Si XII 
brightens significantly and it is common that the intensity of 
the line doubles, and the Ly a fades. This general description is 
based on several studies made on different CME events which 
reached similar conclusion. For instance, broad O VI lines and 
mor e intense Si XII lin e have been observed i n 2000 June 28 
(Ci aravella etaTI l2005h . 20 00 September 12 (ISuleiman et al] 
2005^ and 2000 October 24 (ICiaravella et al.1 12006')^ event. In the 
shock related to the CME event on June 11, 1998 the derived 
compression and ionizati on appear modest wit h respect to the 
inferred shock strength (iRavmond et al.l l2000h . Signatures of 
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shock have been obse rved during the 2000 March 3 CME event 
dMancuso et al.ll2002h . They pointed out that the Oxygen tem- 
perature (To) was much larger than the proton temperature (Tp) 
and that the coronal O VI emission during a shock is well fitted 
by a narrow component which comes from the quiet corona 
and a broad component which comes from the shocked plasma. 
Actually, particles of different species can be heated by different 
amounts in a collision-less shock, since there is not enough 
time for Coulomb collisions to bring the temperatures into 
equilibrium. Observations of shocks in th e solar wind generally 
show that electrons are weakly heated (ISchwartz et al.l l2000l) 
and that minor ions of mass rrii have higher tem peratures than 
proto ns, perhaps as high as = (mi/mp)Tp dKorreck et al.l 
I2007h . Observations of supernova remnants instead show 
Tf, ^ Tj for relatively sl ow shocks and <C Ti for fast shocks 
dGhavamian et al.l 1200 lb and Tp ^ for slow shocks and 
To ~ 8Tp for a fast one dRavmond et alj|2003bl; iKorreck et alj 
12004 . 

Since the whole sample of events amounts to ^ 10 shocks 
linked to CMEs and only few of them have good enough data 
to do detailed analysis, questions about the line emission from 
shocked plasma are still open. Overall it is important to link all 
the observed features listed above with models, at least quali- 
tatively. Because of the large number of parameters that could 
characterize the evolution of a CME it is quite impossible and 
even useless to model quantitatively one single event. On the 
other hand, many recent efforts have focused on the theoretical 
modeling of the solar corona, with the goal to properly describe 
the processes that lead to eruption and activity in the corona. 
We think that modeling should be linked as much as possible 
to observables, because several assumptions lay between the 
physics that the model describes and the actual emission that 
could be inferred from that model. As an immediate implication, 
any attempt to model observables in shock wave propagation 
cannot neglect that the highly dynamic plasma c an be easil y 
in Non Equilibrium of Ionization (NEI) ( Spadaro et al J [19941) . 
since the time scales are too short to allow the ionization state 
of the plasma to relax. 

The aim of this work is to find unambiguous spectral shock sig- 
natures, and to provide a guide for the interpretation of general 
features that could appear in UVCS observations. Our approach 
is to model in detail and diagnose the propagation of a shock 
wave generated from a supersonic fragment of a CME in the 
magnetized coro na. We use a numeric al Magnetohydrodynamic 
(MHD) model dPagano et al.1 l2007b to describe the MHD 
evolution and from the result we compute the plasma emission 
for the O VI and Si XII lines which are diagnostically important 
in UVCS observations of shock waves including the effects 
of NEI. As a guideline to interpret the results, we support the 
study with a simple analytical model of a Rankine-Hugoniot 
planar and adiabatic shock wave. From the radiation model we 
synthesize observable quantities, such as the profiles of spectral 
lines. 

In Section |2] we describe our model, Sec|3] includes the results 
discussed in Sec|4] 

2. The model 

We use a full MHD model of the solar corona tuned to inves- 
tigate the coronal emission visible from UVCS during a shock 
wave propagation. We model a supersonic CME core fragment 
moving upward in a magnetohydrostatic solar corona. During 
its propagation, the fast cloud generates a series of shock waves, 
and we study how the waves perturb the quiet corona. To this 



aim, we need to model the evolution of both the plasma and the 
magnetic field and, to compute emission, also of the ion abun- 
dances. 

We consider the sol ar gravity, the radiative losses (e.g. 
iRavmond & SmithI dl977 l). a phenomenologica l coronal heat ing 
term and the field-oriented thermal conduction C SpitzeJ 19621) as 
important physical effects. The full MHD equation we solve are: 



^4 

dt 

dpv 
du 
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with the constraint given by: 
V • B = 
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where t is the time, p is the density, n the number density, p 
the thermal pressure, T the temperature (T^ — Tp = Ti in 
this model), v the plasma flow speed, u the total energy (inter- 
nal energy E plus kinetic), g the gravity acceleration, Fc is the 
conductive flux according to Spitzer (1962) and corrected for 
the saturation effect (Cowie & McKee 1977), P{T) the radia - 
tive losses per unit emission measure (Raymond & Smit hiil977h . 
Hq = nQP{To) is a constant heating term whose only role is to 
keep steady the unperturbed corona by balancing exactly its ra- 
diative losses, with npfr) = n{r,t = 0) and Tq = T{t — 0) 
(iPagano et al.l 1200 7*). We do not include plasma resistivity ef- 
fects, which can be considered globally negligible on large 
scales in the solar corona. 

We solve numerically the set of the ideal full MHD 
equations with the MHD module of the advanced par- 
allel FLASH code, basicaUy developed by the ASC / 
Alliance Cen ter for Astroph ysical Thermonuclear Flashes in 
Chicago (Frv xell et al.l 20001). vyith A daptive Mesh Refinement 
(PARAMESH, 'MacNeiseeTinilOOO). We include the FLASH 
module for the anisotr opic thermal conduction (Spitzer 1962|) 



implemented by Pagan o et al. ( 20071). and improved for the sat- 



uration effects by Orlando et al. 

We compute the ionization state of the plasma from the his- 
tory and distribution of v, and T obtained with the MHD 
model. This is done asynchronously of the MHD computations; 
we can do this safely because we sample the solution on time 
bins (~ 25 s) shorter than a typical UVCS exposure time (^ 120 
s) and than the shock passage time scales 500 s for a scale 
length of 5 X 10^ cm and a typical sound speed of 10^ cm/s) 
and because we do not expect effects from small scale mix- 
ing by turbulent motion in the presence of the magnetic field 
and of the thermal conduction. The ionization state is computed 
considering the lagrangian transport of ions, and the ioniza- 
tion/recombination processes, i.e. we are not assuming ioniza- 
tion equilibrium during the simulation. The variation of the ion- 
ization fraction of each ion species is governed by the equation: 



dt 



-V-nfv+n,[nf_^,af_^,+nf_,qf_,~nfiaf + qf)](8) 
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where nf is the density of the element Z in the i-th ionization 
level in units of cm~^, is the electron number density and af 
and qf are respectively the recombination and ionization rates 
for the element Z in the i-th ionization level in units of err? j s. 

The ionization state evolution explicitly depends on elec- 
tron and ion densities and the bulk velocity of the plasma, and 
implicitly on the temperature through af and qf . In Eq. ([8]l 
the first term on the right hand side is the advection term, and 
the second is the variation due to ionization and recombina- 
tion processes that take place at the temperature of the plasma 
element. We solve this equation numerically for each relevant 
species and ions, assuming initial ionization equilibrium, and 
treating explicit ly the advection term with a Godunov scheme 
(lGodunovlll959[) and implicitly t he se cond term with ionization 
rates given by 'Cox & Raym ond! (11985 ) and recombination ones 
by Bryans et al. (2006). The evolution of the electron density and 
of the plasma motion throughout the domain are taken from the 
MHD model results (sampled every 25 s). A Godunov scheme 
is applied to compute the ion transfer through the cells due to 
the bulk motion of the plasma, and the ionization fraction is up- 
dated considering the ionizing and recombining collisions (sec- 
ond term in Eq. (O) between electrons and ions. 

The intensity of the lines is computed as the sum of the colli- 
sional (/c) and radiative (Ir) contributions (see the Appendix lAl 
for more details): 

I^h{7ie.nf,T)+Ir{nf,T,v) (9) 

We model the evolution of a dense cloud moving upward su- 
personically in the outer coronal atmosphere as shown in Fig. 
[U The center of the cloud (i.e. the CME fragment) is initially 

Log1U( lemperature} LoglU(Uensity number) 
Time = s 




Fig. 1. Initial condition of the model. The panel shows a section 
at y = of the temperature and density number spatial distribu- 
tions. The the magnetic field lines are also shown (white lines). 

positioned at O.15i?0 above the photosphere and has an upward 
initial velocity. It is denser than the surroundings, but in pressure 
equilibrium with them (i.e. colder than the surroundings atmo- 
sphere). The atmosphere is isothermal at 1.5 MK and stratified 
by gravity in such a way to have n = IQ^cm^^ at the initial 
CME fragment center height. The background magnetic field is 
modeled as a dipole laying in the center of the Sun and seen from 



the central axis. On this basic configuration, we build various se- 
tups of the other parameter and initial conditions. As reference 
set up, the CME fragment has an initial velocity Vc = 1000 km/s 
and is 10 times denser than the surroundings (i.e. ric/riQ — 10 
and temperature: Tc = 1.5 x 10^ K). The intensity of the dipole 
is chosen to give /3 ~ 1 at the initial height of the CME fragment 
center. The other simulations that we present differ only for one 
parameter, and in particular, (3 = \Q, (3 — 0.1, ndriQ — 4, 
Vc = 700 km/s 

Simulations /301 and /310 are devoted to investigate the ef- 
fect of the magnetic field on the global structure of the shocks, 
i.e. whether a stronger (weaker) magnetic field could lead to 
less (more) compression and, ultimately, to less (more) emis- 
sion. Simulation rifA and 1)700 explore situations in which the 
initial shocking fragment carries less momentum (less mass the 
former, less velocity the latter). 

The simulations are all performed in a 3-D cartesian domain 
(x,y,z) that is (8 x 8 x 16) x lO^*' cm large, for the duration 
needed by the shock to propagate well above 2Rq which is the 
height where we put a hypothetical UVCS observations. The 
time is 1000s when the initial velocity of the fragment is 1000 
km/s and ^ 1500s when it is 700 km/s. 

3. Results 

We now discuss our reference simulation and later the others 
which differ from it by just one parameter of the initial condi- 
tions (i.e.: /?, initial velocity of the cloud, density of the cloud). 

3. 1. The reference simulation 

In Fig. |2] the density and the temperature of the plasma after 
1000 s of evolution are shown (for the reference simulation). 
In all the simulations shock fronts depart from the upper part 
of the high speed cloud. Here we address the shock evolution 
only, and not the cloud evolution. The cloud continuously shocks 
the surrounding medium, since its speed remains faster than the 
sound speed (at least during the time of our simulation), it cools 
down during the motion because of the radiative losses and its 
core becomes thermally unstable. 



Log! U( I emperature) LoglU(Uensity number) 
Time = 1000 s 




Fig. 2. As in Fig.[T] but at < = 1000 s in the reference simulation. 
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Table 1. Numerical simulations 



Name 


fJ 


Uc/no 


Vc (km/s) 




1 


10 


1000 


poi 


0.1 


10 


1000 


(310 


10 


10 


1000 


nA 


1 


4 


1000 


v700 


1 


10 


700 



3.1.1. The Shock 

In the reference simulation, the shock propagates radially from 
the cloud. Several shock fronts are generated during the passage 
of the cloud and finally they merge in a bow shock edge ex- 
panding in time. After 1000 s of evolution the cloud has reached 
the height of ~ 2.5Rq and the produced shock edge has trav- 
eled ^ 0.4i?Q from the originating cloud. Behind the edge of 
the shock wave, the structure of the shocked gas is due to the 
interaction and superposition of the different shock waves. 

Since this work addresses the signatures of a shock transi- 
tion during a CME, we focus our attention on the shocking edge 
which drives the disturbance of the quiet outer corona rather than 
on the inner structure of the shock. The closer in x direction the 
shocked plasma is to the shocking cloud, the less the shock has 
propagated in the corona, the stronger is the shock, and the more 
the bulk velocity of the post-shock plasma is perpendicular to the 
solar surface. The cloud velocity remains superalfvenic and the 
magnetic field is shocked as well. For this reason the post-shock 
magnetic field tends to be parallel to the shock front. The plasma 
P is relatively high in the upper corona 15 in our simulation 
at 2Rq), and so the magnetic field pressure is negligible with re- 
spect to the plasma compression, while the magnetic field orien- 
tation influences the thermal properties of the shock and, as con- 
sequence, the ionization state. Since the shocked magnetic field 
is parallel to the shock front which is hotter than the surrounding 
atmosphere, the thermal conduction toward pre-shock regions 
is ineffective and the shocked plasma does not cool down. The 
shock heats the plasma, the ionization is enhanced and the ion- 
ization state changes. Immediately behind the shock the plasma 
is far from the ionization equilibrium at the shock temperature, 
and the ionization equilibrium is approached farther from the 
shock front, where the temperature has been no longer changing 
much and the plasma is slowly relaxing. 

3.1 .2. O VI and Si XII shock emission 

We apply our general radiation model to O VI and Si XII Unes, 
that are important for the CME shock front diagnostics. Fig. [3] 
shows the evolution of the line intensity across the UVCS slit 
computed from the full MHD simulation and considering non- 
equilibrium of ionization. The intensity jump across the shock is 
sharper in the Si XII line than in the O VI line. Our model shows 
that the shocked plasma is far from ionization equilibrium; in 
fact, the heating time due to the shock is much shorter than the 
ionization equilibrium time and there is not enough time for the 
ionization state to change as soon as the shock has reached the 
plasma. For instance, the ionization time for the O VI at the post- 
shock temperature (T ^ 3 MK) is roughly Teg ^ 400 s (see 
AppendixlB), significantly longer than the heating time. We can 
estimate the actual heating time {Theat) from the speed of the 
shock and the gyroradius of the particles. Considering a mag- 
netic field of B 0.5 G, a temperature of ~ 2 MK and the ve- 
locity of the shock 500fcm/s) we get Tteat ~ 10~^ s. Since 
the temperature in the region behind the shock is higher than the 
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Fig. 3. Evolution of the emission of O VI (upper panel) and Si 
XII (lower panel) expected along an UVCS slit positioned at 
2i?Q integrated along the line of sight. The evolution is sampled 
every 25 s. The white lines mark the shock edge. 



peak temperature for the O VI, the emission becomes fainter and 
fainter behind the shock while the ions approach equilibrium, i.e. 
at later time. The shock heating makes the Si XII abundance in- 
crease, because the Si XII peak temperature is higher than the 
temperature of the quiet atmosphere (1.5 MK). Because of this, 
a large fraction of the Si XII emission comes from the inner part 
of the shocked region. 

As mentioned above, for diagnostic purposes, we focus the 
analysis at a characte ristic height for the UV CS observations and 
for shock formation (iRavmond et al.ll2"003al) . i.e. 2i?0. It should 
be noted, that, in spite of the formation of the shock, the CME 
may either accelerate or decelerate that far from the solar sur- 
face and that long after the ignition phase. In order to investigate 
the role of the non-equilibrium ionization state, of the thermal 
conduction, and of the magnetic field, we compare the intensity 
of the plasma behind the shock obtained from the MHD simula- 
tion with the one predicted by a simple model of adiabatic and 
locally planar shock according to the Rankine-Hugoniot jump 
relations of density, temperature and velocity as a function of 
the post-shock velocity of the plasma. For further details about 
this analytical model see the Appendix |B] This simple model 
is a rough approximation for our shocks, since the shock is not 
propagating in a very low beta plasma (neglect the momentum 
flux due to the Lorentz force), the thermal conduction is mostly 
inhibited by the magnetic field (adiabatic plasma) and the shock 
heating is much faster than the plasma ionization and recombina- 
tion (the ionization fraction does not change at the shock front). 
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Significant departures from this approximation may be signature 
that some of the assumptions listed above do not hold. 

To compare the detailed simulation with the analytical model 
so as to evaluate the deviations, in Fig.|4]we plot the ratios be- 
tween the post-shocked and pre-shocked density and tempera- 
ture as function of the post-shock velocity of the plasma for the 
MHD simulation and for the analytical model. The post-shock 
quantities for the simulation (i.e. post-shock velocity, density 
and temperature jump) are computed at the shock edge, marked 
by white lines in Fig. [3] 
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n l 

12 3 4 5 

Post-shock plasma velocity (10' cm/s) 

Fig. 4. Ratio between post-shock and pre-shock temperature 
(solid lines) and density (dashed) as a function of the post- 
shock plasma velocity. Thick lines mark results obtained from 
the MHD simulation, and thin lines from an analytical model of 
Rankine-Hugoniot relations 

on an adiabatic and planar shock. 

The front of the bow shock begins to cross the UVCS slit 
field of view at ^ 675 s. UVCS detects first the upper and 
fastest part of this front and later the flanks. As mentioned above, 
the shock strength decreases in time during the crossing of the 
UVCS slit field of view and at t=700 s the plasma is acceler- 
ated to ~ 500 km/s because of the shock and to ^ 250 km/s at 
t = 1000 s (see Fig.|4]i. As shown in Fig.|4]the density jump in 
the simulation roughly agrees with that of the analytical model, 
meaning that the magnetic field has little relevance in this dy- 
namic regime. The temperatures agree with less accuracy be- 
cause the thermal conduction is not completely ineffective, and 
the shock can cool down for that. For instance, the thermal con- 
duction is effective in the region where the shock is stronger, i.e. 
where the shock propagation is not perpendicular to the post- 
shock magnetic field. Here (i.e. x ^ 2 x 10^" s and z ^ 10x10^" 
cm in Fig. |2]i, the angle between the magnetic field and the ther- 
mal gradient is ^ 30°. Considering the thermal gradient, the 
temperature and the length scale involved, we estimate a thermal 
conduction time scale of ^ 5 s to be compared to the dynamical 
(expansion) time on the same length scale which is ^ 100s. The 
thermal conduction is therefore very effective and leads to tem- 
peratures lower than those predicted by the Rankine-Hugoniot 
relations. 

From the jump in the density and the temperature we derive 
the jump in the emission for O VI and Si XII with our radia- 
tion model. It is plotted in Fig.|5]as a function of the post-shock 
plasma velocity (O VI upper panel and Si XII lower panel). In 
the O VI line the radiative part of the line is significant for the 
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Post-shock plasma vslodty (1 0' cm/s) 

Fig. 5. Ratio between post-shock and pre-shock emission in the 
O VI (upper panel) and Si XII (lower panel) lines as a function 
of the post-shock plasma velocity. Results from the MHD simu- 
lation (thick lines) and from the analytical model (thin lines) are 
shown. In the latter case we consider approximated ionization 
state formulae assuming r = (solid line) or r = 100 (dashed 
line) and our radiation model (see AppendixlAb. 



unperturbed plasma, but it completely vanishes because of the 
Doppler dimming in the region behind the shock. The radiative 
contribution to O VI 1032 is 4 times that of O VI 1037 (in the 
absence of other pumping), while the collisional component is 
a factor of 2 brighter. The shock compression makes the colli- 
sional contribution of the line increase, whereas the shock ac- 
celeration makes the radiative contribution decrease at the same 
time. For this reason, the shock shows an emission fainter than, 
or at m ost comparable with, the qu iet corona, as usually ob- 
served (ICiaraveUa et alJl200l l2006h . i.e. O VI emission ratio 
'--^ 1 in Fig. |5] The weight of the radiative part depends on the 
non-shocked plasma density and in our model considering an 
unperturbed density of ~ 2 x lO^cm^'^, it is 85% of the emis- 
sion and increases as the unperturbed density decreases. Instead, 
since the radiative part of the Si XII is negligible, the total in- 
tensity of the post-shock line emission jumps significantly and 
relatively to the strength of the shock, as plotted in Fig. |5] (lower 
panel). As a consequence, the Si XII line clearly b rightens, while 
the O VI line does not, as lCiaravella et al.l (12006 ) pointed out. 
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3.2. Comparison with otiier simulations 
3.2.1. Different (3 

As anticipated above, we perform two more simulations con- 
sidering different intensities of the magnetic field, with a less 
(/3 ~ 10) and more (/3 ~ 0.1) intense magnetic field than in 
the reference simulation, respectively. Fig. |6] shows the density 




Post-shock plasma velocity (1 0' cm/s) 



Fig. 6. Ratio between post-shock and pre-shock temperature 
(solid lines) and density (dashed lines) as function of the post- 
shock plasma velocity. The thickness of the line is inversely pro- 
portional to the initial /3 of the simulation, (P — 0.1, (3 = 1 and 
P = 10). The thinnest lines is for the ratio computed analytically 
from Rankine-Hugoniot relations 



and temperature jump across the shock for the simulations with 
different values of (3. The simulation with /? = 10 shows a very 
similar behaviour and the shock properties do not depart signif- 
icantly from the reference case already discussed. This is not 
surprising, since we have already shown that the magnetic field 
has no relevance for the shock dynamics in this regime. 

In a stronger magnetic field (Fig.|6]l, the shock is significantly 
influenced by the magnetic pressure. According to the theory of 
MHD shocks the compression due to the shock is splitted be- 
tween magnetic field and plasma and, as result, the plasma den- 
sity jump is weaker than in a normal HD shock ( Baz er & EricsonI 
Il959i) . This effect characterizes our simulation with a stronger 
field. The density jump from the pre-shock to the post-shock re- 
gion is 5% less than reference case, but the temperature does 
not show any significant departure since the energetics of the 
shock is not influenced by the magnetic field except for the par- 
tial thermal insulation which acts regardless of the f3 value. The 
radiative properties of this case descend from these considera- 
tions. The emission of the O VI and Si XII line is weaker than 
in the reference case, because the density is lower. The radiative 
part of the line vanishes because of the Doppler dimming, and 
the intensity of the line is weaker than the reference case is by 
- 10%. 



3.2.2. Weaker cloud 

We now discuss the propagation of shocking clouds that initially 
carries less momentum, i.e. with initial density contrast n/nO = 
4, instead of n/nO = 10, or initial cloud velocity vq = 700 km/s, 
instead of vq ~ 1000 km/s. The shock propagation is slower in 



those simulations and the top of the cloud needs 600 s and 875 s 
to reach the height of 2Rq, in the simulation with n/nO = 4 and 
in the one with vq = 700 km/s respectively. During the travel 
of the cloud in the corona, the shock departs from it radially 
reaching a final distance from the cloud of ~ O.6i?0 in the x 
direction. The final bow shape of the shock front is similar to the 
one in the reference case. Of course, the smaller the energy of the 
cloud, the less the shock strength at any time, so that the plasma 
is never accelerated over ^ 300 km/s. The velocity relative to 
the cloud is 30% smaller and that leads to a corresponding delay 
in the detection at 2Rq. This delay allows the plasma to better 
approach ionization equilibrium. 

3.3. O VI line profile 

To further investigate the observational constraints for the shock 
properties we show here the line profile of the O VI lines synthe- 
siz ed from our reference model. This line is analyzed in detail 
by iMancuso et aL I (I2002h who indicate shock diagnostics from 
UVCS observations. Here we try to provide further insight to 
that analysis, but our synthesis of the line profile is very gen- 
eral and could be applied also to the Si XII line after tuning the 
appropriate parameters. 



_Pre-shock (t = 500 s) Line profile 
_Post-shock (t=900 s) Line profile 

from M>1 

from M<1 




lOJO lOJl 1032 ^CSS 10J4 

Angstrom 

Fig. 7. Line profile of the O VI line synthesized from our ref- 
erence model if observed as a limb CME at t ^ 500 s (point- 
dashed thick line) and at < 900 s (solid thick line), i.e. before 
and after the shock crossed the UVCS field of view, respectively. 
The emission is averaged on time bins of 150 s and on space 
bins of 9' along the UVCS slit through which the CME is ob- 
served. The profile is normalized to the peak. We also show the 
two components which combine to form the line at t ~ 900: the 
one emitted by the shocked plasma {A4 > 1, solid thin line) and 
the one emitted by the unshocked plasma (Ai < 1, dashed thin 
line). 



We take into account the thermal broadening of the line and 
the Doppler shift due to the plasma motion projected along the 
line of sight, but not the instrumental broadening. To approach a 
realistic UVCS observation we average the spectrum over time 
bins of 150 s and over space bins of 9'. We also consider in- 
tegrations along two different lines of sight: one mimicking a 
limb CME, the other a halo CME. Finally we discuss the effect 
of the differential Oxygen/protons shock heating on the Limb 
CME line profile. 
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3.3.1. Limb CME 

In Fig. |7] we plot the profile of the O VI Une for the reference 
model as if it were observed as a limb CME before and after the 
shock passed the UVCS slit. For the moment we assume equal 
proton, electron and ion temperature. After the shock the central 
part of the line has not changed because it is emitted by that part 
of the plasma never involved in the shock passage. Evidence of 
the passage of the shock across the UVCS slit can instead be 
found in the wings of the later line. They are more prominent 
than in the quiet line, because of the line shift due to the line of 
sight component of the shock velocity. The resulting line profile 
is the sum of the emission contribution from un-shocked plasma 
and shocked plasma. The un-shocked plasma emits only a nar- 
row component (Fig.|7] Ai < 1 line), while superposition of the 
lines emitted by the shocked plasma gives both a central com- 
ponent and significant wings. These two emission contributions 
merge in a two component line. The first is a naiTow and bright 
component, the sharp peak emitted by the unshocked corona and 
by the shocked plasma that moves perpendicularly to the line of 
sight; the second component is due to the shocked plasma that 
has a significant component of the velocity parallel to the line of 
sight, so that it emits thermally broadened lines centered far from 
1032A. The relative weight of the two components is strongly 
sensitive to the backgr ound emission, which seriously influence 
the narrow component. i Mancuso et al.l ( |2002|) indeed observed a 
similar line profile and concluded that a shock was crossing the 
field of view of the UVCS slit during the observation. We argue 
that a two component line profile, as shown in Fig. |7] may in 
general be a signature of the presence of a shock. 

3.3.2. Oxygen shock heating 

We expect that the wings of the line profile become more promi- 
nent if the shock heating mechanism is more effective on the 
heavy ions than on the protons, because the thermal broad- 
ening becomes larger right in the Doppler-shifted components 
of the ion lin e s. To test this consideration, already made by 
iKorreck et all (|2004|) and [Raymond etal] (|2003b), we check 
what happens if we artificially increase the O VI temperature 
obtained from our MHD model by a given factor. Fig. |8] shows 
the O VI li rie profile, obtained by increasing the tem perature by 
a factor 8 ('Korre ck et ani2004t iRavmond et al.l | 200 3b) and 16 
(oxygen atomic mass). Because of the temperature decoupling 
between electrons and ions, the FWHM of the line emitted by 
the shocked plasma increases from ~ 0.2A up to ~ 1 A. 

The wings of the line profile brighten if a greater Tqvi /Tp 
ratio is supposed. Fig.|9]shows the ratio between the total emis- 
sion of the wings (i.e. emission out of the range 1032 ± O.SA) 
and that of the center of the line (i.e. inside that range). It in- 
creases from 0.1 {To vi = T^) up to 0.24 (Tq y/ = 16Te). The 
trend is mostly due to the increasing thermal broadening of the 
emission from the post-shock plasma (Fig. |9] dotted line), but 
it is limited by the presence of the other broadening component 
which is independent of the temperature (Fig.|9] thin solid line). 
The sum of the two effects gives the total strength of the wings 
in the observed lines (Fig.|9] thick solid line). The stronger is the 
shock, the stronger are the wings. Therefore at ~ 575 s, when the 
first part of the shock crosses the UVCS slit field of view and the 
shock is the strongest, the wing/center ratio is ^ 0.6 and later, at 
~ 1000 s, it decreases down to 0.24 (shown in Fig.|9j. It should 
be noted that Fig. [8] shows that Tn/Tj, = 16 resembles the 
profile s shown by [Raymond et al.l (l2003bl) and iMancuso et al.l 
(I2OO2I) far better than the lower temperature ratios, supporting 
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Fig. 8. As in Fig.|7] but considering different Oxygen-electron 
plasma temperature ratio, as indicated in the legend. All the line 
profiles are computed att ^ 900s. 
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Fig. 9. Fraction of the O VI line emitted in the wings as function 
of Oxygen-electron plasma temperature ratio. The thin solid line 
is considering only the Doppler effect (temperature of the atmo- 
sphere not changed), the dashed line only the thermal broaden- 
ing, the thick sohd fine both of them. 



the assertions of those papers that Tp/Tp was large. Our ap- 
proach, when applied to specific events, can give reliable and 
alternative measurements on the actual Ti/Tp, once measured 
the post-shock electron temperature, and of the components of 
the velocity, since the value of the wing/center ratio uniquely 
corresponds to a certain ion temperature. 

3.3.3. Halo CME 

Fig. [To] shows the results for the reference model if oriented as 
a Halo CME expelled toward the observer, not considering any 
temperature increase for the ions. In the line profile there is only 
one tail on the blue side. The line of sight is appropriate to detect 
the Doppler shift due to the global front motion, rather than that 
due to front expansion, detected for the Limb CME line of sight. 
At increasing distance from the shock central axis the blue-shift 
is progressively reduced due to the projection effect. The total 
line profile consists of the superposition of radiation with con- 
tinuously varying blue-shift. Fig. [TOl shows the line detected by 
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Fig. 10. Line profile of the O VI line synthesized from our refer- 
ence model if observed as an Halo CME at t ~ 900 s, i.e. after 
the shock crossed the UVCS field of view. The emission is av- 
eraged over time bins of 150 s and over space bins of 9' along 
the UVCS slit through which the CME is observed. The profile 
is normalized to the peak. For the solid line the slit is ~ 5' apart 
from the center of the Sun, ^ 7' for the dashed one. 



the slit pointing at ~ 5' (solid line) and ~ 7' (dashed line) from 
the center of the Sun. The closer the field of view to the Sun 
disk, the larger is the line of sight component of the post-shock 
bulk velocity and the more blue-shifted is the emission. The line 
presents a very bright quiet emission and an envelope of shifted 
lines extending to ~ lA from the quite line (see Fig.[TOl solid 
line). Although the details of the structure of the line are strongly 
related to the structure and geometry of the shock, the very pres- 
ence of this smooth and asymmetric component of the line could 
be a distinctive feature, although such observations are techni- 
cally unfeasible for UVCS, because the slit is too close to the 
center of the Sun, we may equally hope to detect this feature in 
real observations. In fact, since the shock front is more than 16' 
wide, if the CME is not expelled perfectly along the Sun-Earth 
line, part of the shock front could be visible out of the Sun disk 
in a realistic UVCS field of view. This conclusion would apply 
to the halo CME shocks observed by ICiaravella et alj (l2()06b . 
Moreover, as we akeady claim above (Section [3.2.21 ). a shock 
front could extend over a wider angle if the shocking fragment 
is less dense. 

Finally, for CMEs expelled in directions in between the ex- 
treme cases considered here we expect to detect both the broad- 
ening effect due to the thermal broadening of the shocked plasma 
(Limb CME) and the one due to the Doppler shifted emission 
from fast plasma moving toward the observer. 

4. Discussion and conciusions 

In this work we study the shock fronts departing from super- 
sonic CME fragments with detailed MHD modeling. We address 
the diagnostics obtainable with UVCS observations, which have 
been widely analyzed in this framework, and in particular the O 
VI and Si XII emission lines. As a result, we are able to explain 
most of the spectral signatures of the shocks revealed by the ob- 
servations, and we step forward to predict more shock properties. 
In particular we indicate the main differences in the line shapes 
expected between halo CMEs and limb CMEs and the strength 



of the wings of the O VI line as possible diagnostics for ion 
shock heating. 

We have chosen as a reference case the shock front originat- 
ing from a CME fragment 10 times denser than the surrounding 
1.5 MK magnetized (/? w 1) corona expelled at a speed of 1000 
km/s. The ambient magnetic field is a large scale dipole. Then 
we have explored other cases with higher and lower magnetic 
field, with less dense cloud and with slower cloud. For each sim- 
ulation we have computed the emission in the selected lines and 
derived the line profiles with a detailed radiation model includ- 
ing the radiative and collisional contributions and the effects of 
non-equilibrium of ionization. 

From the reference simulation we find that some of the dis- 
tinctive shock features detectable in O VI and Si XII lines arise 
from the relative brightness of the collisional and radiative con- 
tributions in the unperturbed corona. Since, the radiative con- 
tribution is actually negligible for Si XII, but important for O 
VI, the intensity jump across the shock is sharper in the Si XII 
line than in the O VI line. The intensity of the Si XII line scales 
roughly as the square of the compression, while the one of the O 
VI scales more smoothly because its emission is strongly biased 
by the radiative contribution. In the O VI line the radiative part 
completely vanishes in the region behind the shock because of 
the Doppler dimming, while the shock compression makes the 
collisional contribution of the line increase. The two effects are 
competitive and comparable, so that the shock shows an emis- 
sion fainter than or comparable with the quiet corona in this line. 
In the Si XII line, the shocks produce only the growth of the col- 
lisional contribution. This explains why at the same time the Si 
XII line clearly brightens, while the shock has little influence on 
the O VI line. These c onsiderations explain past observations. 
ICiaravella et alj (l2005h found that during the shock passage in 
the 2000 June 28 event the O VI em ission was weaker or c om- 
parable to the quiet corona emission. ICiaravella et al.l(l2006l) ex- 
amined several events and showed that it occurs commonly that 
when a shock is detected the O VI does not brighten, and the Si 
XII brightens significantly. 

The ionization fractions of oxygen and silicon do not change 
in shock transitions because the shock heating and passage are 
much faster that the ionization/recombination processes. In our 
simulations the shocks are not extremely strong, so that they 
represent realistic situations, typically observed. For this rea- 
son, we expect that the conditions of ionization non-equilibrium 
are commonly found and that the ionization fractions are com- 
parable to those of the quiet corona. Behind the shock front 
the plasma should slowly approach the ionization equilibrium 
at post-shock temperature. 

The effect of the magnetic field on the shocks depends on 
the strength of the field. In high f3 plasma the role of the mag- 
netic field is only to thermally insulate the shocked plasma. It 
happens because in this regime the shock is superalfvenic and it 
is responsible for the reorientation of the magnetic field parallel 
to the shock front, i.e. perpendicular to the thermal gradient. In 
low /3 plasma, the magnetic field influences also the compression 
of the plasma. In this regime (also superalfvenic) the compres- 
sion due to the shock is split between the magnetic field and the 
plasma; therefore, the de nsity (and the li ne intensity) is lower 
than in a high f3 regime. iRavmond et alj ('2000) argued that a 
relatively strong magnetic field could in principle be responsi- 
ble for modest compressions with respect to the observed shock 
speed. 

Differences in the kinetic energy of the CME could lead to 
different extent of the shock front. We find that the weaker the 
CME core, the more extended and slow is the shock propagation. 
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When we observe the shock propagating from a limb CME 
we expect to observe a profile of the O VI line consisting on two 
distinct overlapping contributions. The unperturbed plasma that 
lies along the line of sight emits a narrow line and the shocked 
plasma e mits a wider line . This combinat ion is the kind of line 
shape that lRavmond et alj (^000) and Ma ncuso et all (|2002|) ob- 
served during the CME events. They used the double component 
line profile as signature of the shock. In our simulation we get 
the same line profile shape after the shock passage. Therefore 
we confirm that this spectral feature is a reliable signature of the 
shock presence. 

In our simulation we have a line of sight component of the 
velocity of ~ 200/cm/s, leading to a line shift of 0.7 A. 
The shock heating Tq — Tp leads to a line shift of ~ 0.2 A. 
The acceleration of the plasma is responsible for the presence of 
the wings and higher Tp/Tp leads to more significant wings. 
For this reason, diagnostic of the ion shock heating could be 
based on measuring the relative importance of the wings of the 
O VI line and on estimating the post-shock electron tempera- 
ture, and the components of the velocity. Here we present an 
example of proton and oxyge n enhanced heating jKorreck et alj 
120041: [Raymond et al.ll2003bl) . It will be interesting in the future 
to apply this kind of diagnostics of the ions shock heating to one 
particular event. 

Finally, we expect to see an asymmetric O VI line profile 
with a smooth tail on the blue-side, when the shock propagat- 
ing from a halo CME expelled toward the observer enters in the 
UVCS field of view. A condition for the detectability of this ef- 
fect out of the solar disk are a large line of sight component of 
the plasma velocity and a wide shock front. 
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Appendix A: The radiative and coiiisionai 
contribution 

The intensity of the lines is modeled as the sum of the coiiisionai 
(/c) and radiative (1^) contributions as functions of iie, nf , T 
and v : 



I^Ieine,nf,T) + Ir{nf,T,v) 

A. 1. Coiiisionai contribution 
The coiiisionai contribution is: 

Icine,nf,T) = UenfQiT) 



(A.l) 



(A.2) 



where Q{T) is the coiiisionai excitation rate coefficient for 
the ions which treats the coUisionally induced emission of the 
ions. Here we present the approximated formula for Q(T) we 



use for the two lines of interest for this paper (i.e. O VI 1032A 
and Si XII 499A): 



QiT) 



O VI 



o y/' 



' /T 



In 



-O VI 



QiT) 



Si XII 



^Si XII ^ 



- 



T 



■In- 



khT 



XII 



(A.3) 



(A.4) 



where Q'^ and Q^"^ are two normalization factors 
useful to match the emission at the ionizatio n equilibrium wit h 
the one prescribed by the CHIANTI database (iLandi et al.l2006l) . 
T^t^-f = 139.000 K and T^Jt^" = 288.000 K are the temper- 
ature for O VI and Si XII respectively and Eq vi = 11. 99 eV 
and £si xii = 14.8 eV are the threshold energy respectively for 
O VI and Si XII. 



A.2. Radiative contribution 



The radiative contribution is: 

Ir{nf,T,v)=ro^faV 



(A.5) 



where /q is the intensity emitted from the Sun disk that 
reaches the height r to be scattered by the coronal ions, a is 
the cross section of the scattering, and V is the Doppler dim- 
ming factor. The radiative contribution is computed only for O 
VI, since it is negligible for Si XII. /q is computed by: 



jr 



tO VI 1032r,^ 

Iq 27r 



1 



1 



-n,Q 



(A.6) 



where i"32 



1.94 p/iot/ (cTO^ s sr) is the disk inten- 
sity and the following function of r is simply the dilution factor 
due to the distance from the disk. 
a is computed (for Oxygen) by: 

o VI 

(A.7) 



VI 

Co VI — CTq 



"disk 



mo 



Where v^J.^ = 30 
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X 10^ cm/s is the wideness of the 
Oxygen lines emitted from the disk, and mo — 2.67 x 10^^"^ 
is the Oxygen mass. 

T> is computed (for Oxygen) by: 

L,0 VI > . ^Lyl3 ^ > 



V. 



O VI 



O VI 1032, 



O VI 1032 



(A.8) 



where the second term accounts for pumping of O VI by 



4.13 X 10 
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Ly/3 (iRavmond & Ciaravellal l2004l) and /g 
phot/(cm^ s sr) is the disk intensity for the Lyman /3 line, 
Vq^^ ~ 18 X 10^ cm/s is the distance between the center of 
the O VI 1032 A line and the LyP line, Vrad is the radial veloc- 
ity of the plasma and Vcut is the cutoff velocity for the Doppler 
dimming which is given by the overlapping between the disk line 
profile and the coronal hne profile. For Oxygen and Hydrogen it 
is respectively: 



V. 



OVI _ , In VI 



cut 



disk 



3hT 
mo 



disk 



3hT 
rriH 



(A.9) 



(A. 10) 



where v^^^^ = 112 x 10^ cm/s is the width for Hydrogen 



lines emitted from the disk, and mn 
Hydrogen mass. 



1.68 X 10"^" g is the 
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Appendix B: Analytical shock model 

From the Rankine-Hugoniot conditions for a planar adiabatic 
shock and our radiation model presented in Appendix|A] we de- 
velop an analytical model for the jump in the intensity of the 
emission lines due to the shock. 

First, we compute the jump in He, nf , T and v as function 
of the Mach Number, then we use the radiation model presented 
in AppendixlAlto get the jump in the intensity of the emission. 

Hereafter we indicate the pre-shock and post-shock quanti- 
ties respectively with the subscript and 1. The jump in elec- 
tron density, temperature and velocity due to the shock are 
jLandau & Lifshitzll 19661) : 



riei 



(7 + l)X2 



neo 2+(7~l)A^2 



Tn 



vi = 2csQ 



(7 +1)2^2 



(7 + l)M 



(B.l) 



(B.2) 



(B.3) 



where A4 is the Mach number and Cg is the sound speed. 

The ionization fraction is computed from the following ap- 
proximated formula (Eq. ( IB.4I )) in which pre-shock ionization 
equilibrium is assumed. 



ff = Ff{Ti) + {Ffin) - Ff{T,))e- 



(B.4) 



where ff is the ionization fraction for the ions Z in the i-th 
ionization state, Ff (T) is the equilibrium ionization fraction for 
the ions Z in the i-th ionization state at temperature T, r is an 
estimation of the time elapsed from the shock front passage, and 
Teq is the time scale needed to get ionization equilibrium which 



is given by: 

Teq = {neQeffY 



(B.5) 



Because a single rate dominates for both of the ions we are 
interested in, we can approximate q^j / by the fastest rate involv- 
ing the ion. 

We can compute the ions density by: 



ft 



(B.6) 



where Af "" is the solar abundance of the ion Z. 

Through this computation we get the values for ne,nf ,T and 
V for the post-shock region, while we know all of them for the 
pre-shock region. At this point we compute the jump in the emis- 
sion lines from the radiation model presented in Appendix lAl 
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